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INTRODUCTION 


As the era of the Space Transportation System progresses toward the possible 
development of a permanent manned presence in space, long-term exposure of astronauts 
and equipment to large fluences of galactic heavy ions indicates a need to investi- 
gate methods of shielding from these high-energy ions* To properly evaluate the 
passive shield requirements, a comprehensive theory describing the interaction and 
propagation of these particles, and their subsequent reaction products, in the space- 
craft structure and in the inhabitants is required. The nuclear interaction model 
used is described in detail in references 1 and 2 (and references therein). The 
heavy-ion propagation is described by an energy-dependent extension of the transport 
theory given in reference 3. Polyethylene and Du Pont Kapton were chosen for the 
analyses since they may be used as constituents of a multilayered shield design to 
replace the usual bulk metallic spacecraft shielding. In the present report, trans- 
port coefficients for arbitrary ions incident upon Kapton and polyethylene absorb- 
ers are presented over a broad range of energies. The results of the transport 
theory are illustrated by calculating the doses from 20 Ne beams at 350, 670, and 
2000 MeV/amu as a function of depth in the absorbers. 


DEPTH-DOSE EXPRESSIONS 

With the straight-ahead approximation and the target secondary fragments 
neglected, the transport equation is written as 



where (x,E) is the flux of ions of type j with atomic mass A^ at x moving 
along the X-axis at energy E in units of MeV/amu, Cj(E) is the corresponding 
macroscopic nuclear absorption cross section, Sj (E) is the specific stopping power 
(change in E per unit distance), and mj^E) is the fragmentation parameter of 
ion j produced in collision by ion k. (A list of symbols appears after the refer- 
ences in this report.) In terms of the specific stopping power, the range of the ion 
is 
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Using an interative procedure (ref. 3), equation (1) is solved by the method of char- 
acteristics. The resultant series solution is used to evaluate the dose as a func- 
tion of depth as follows: 
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for a monoenergetic beam on type-M ions of energy E Q , where the stopping power is 


Sj - A jSj 


The dose from the radiation field of the primary (incident) beam, obtained by setting 
the right-hand side of equation (1) equal to zero, is given by 
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where the nuclear attenuation factors are 


P m (E) = ex P [-0 M (E) R m (E)] 
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The average extinction coefficient is defined as 
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and the residual energy is 
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The first perturbation to the homogeneous solution yields the dose contribution from 
the secondary ions produced by the fragmentation of the primary beam. This contribu- 
tion is 



where the energy spanned by these secondary ions is given by the "lower limit" 
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and "upper limit 
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The range scale parameter for the type-j ion is (ref. 3) 
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The second perturbation, which yields the dose contribution from the tertiary ions, 
is 


d (2 ’(,0 


A.V./e' - E! \ 

u j l 


jk ^kV'kM 0 *! (V M - v k )( a M “ ^ k ) x 


-a. x 
, 0 


- e 


-v 


-a.x 
> J 


- e 


V 


°M 


a. 

3 


o. - a. 
k 3 


(13) 


where the energy range spanned by these tertiary ions is given by the upper limit 
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and the corresponding lower limit 



(15) 


The parameters m and a of equation (13) are evaluated at E Q . The results of 
equations (14) and (15) are understood to be zero whenever the right-hand sides are 
negative. The above expressions can be applied to various shield materials of uni- 
form composition. Each specific application requires knowledge of the appropriate 

transport coefficients S.(E), a., and m. v . 

J J J K 


TRANSPORT COEFFICIENTS 
Stopping Power 

The total stopping power is obtained by adding the electronic and nuclear con- 
tributions according to the detailed methods given in reference 3. The results for 
various ions in polyethylene, as a function of energy, are shown in figure 1 . Fig- 
ure 2 shows the comparable results for Kapton. 


3 



Nuclear-Absorption Cross Section 


The macroscopic nuclear-absorption cross section 
expression 


(E) 


is determined from the 


a. (E) = Tp. a. . (E) (16) 

3 i 1 13 

where the values of are the elemental constituent number densities for the 

absorber , and the values of (E) are the ’Microscopic" nuclear-absorption cross 

sections. Values of are listed in table 1. The values of a ^ were taken from 

references 1 and 2. Table 2 is a listing of the resultant macroscopic cross sections 
a. for various projectile ions, as a function of incident energy, in Kapton. The 
results for polyethylene are shown table 3. The average extinction coefficients for 
several projectiles are shown in table 4 for Kapton and table 5 for polyethylene. 
Values for intermediate energies and mass numbers can be found by numerical 
interpolation. 


Nuclear Fragmentation Parameters 

The basic fragmentation parameters, or multiplicities, for projectile nuclei 
colliding with an absorber are determined from the expression 
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where a (Z^,A^, A T , Z . ,A • , E) is the partial production cross section for a fragment of 
type A. and Z> produced by an ion of type and Z^ colliding with a target 

of mass A t , and where Cj c (A t ,E) is the macroscopic absorption cross section for the 
kth incident ion colliding with the target of mass number A^, • The partial produc- 
tion cross sections are those obtained from the semiempirical formulas of Silberberg 
and Tsao (ref. 4) augmented by the light fragment production cross sections of 
Bertini (ref. 5). These have been modified, however, in two ways: (1) Rather than 

scaling by total ion kinetic energy, as suggested in reference 4, the relative target 
velocity (energy per nucleon) was assumed to be the appropriate parameter for evalu- 
ating the hydrogen fragmentation cross sections, since particle velocity, rather 
than total kinetic energy, is approximately conserved in these interactions (refs. 6 
and 7); and (2) To account for the lack of fragment mass and charge conservation in 
reference 4, the fragmentation parameters are renormalized to ensure mass and charge 
conservation. The multiplicative renormalization factor is 

F = [ (Z A )/(Z A o )] 1/2 (18) 

L p p s s J 

where Z g and A g are the total fragment charge and mass obtained from the formula- 
tion of reference 4, and Z and A are the incident projectile ion charge and 
mass. Fragmentation parameters and total renormalized fragment charge for 20 Ne beams 
onto Kapton and polyethylene are given in tables 6 and 7 for a broad range of 
energies • 
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Ne DEPTH-DOSE RELATIONS 


20 


90 

The relative doses as a function of depth for Ne beams at 350, 670, and 
2000 MeV/amu, obtained from equations (5), (9), and (13), are displayed in fig- 
ures 3 through 5 for a Kapton absorber, and in figures 6 through 8 for a polyethylene 
absorber. The curves at 350 and 670 MeV/amu show that the dose from secondary and 
tertiary ions is substantial (nearly 50 percent of the entering dose) for several 
centimeters beyond the depths where the primary beam is completely attenuated (in 
the region of the Bragg peak). At 2000 MeV/amu no Bragg peak is visible. The large 
production of secondaries and tertiaries, however, yields substantial doses even 
though the primary beam is largely attenuated. For a Kapton absorber, for example, 
figure 5 indicates that 30 cm of absorber reduces the dose from the primary beam to 
about 10 percent of its initial value at the absorber surface. The production of 
secondaries and tertiaries due to nuclear fragmentations, however, increases the 
total dose to over 25 percent of its initial value. It is not until an absorber 
depth of about 50 cm is reached that the total dose is reduced to 10 percent of its 
initial value. At this depth the primary beam dose is only 2 percent of its initial 
value. Comparable results are seen in figure 8 for polyethylene. A relative com- 
parison of the effectiveness of Kapton and polyethylene as a shield material is not 
included in this report, because these predicted depth -dose values have not yet been 
validated by experiment. In addition, a more meaningful comparison needs to consider 
the isotopic composition of the total beam, since different nuclear species produced 
by the fragmentations could have drastically different biological effects on astro- 
nauts because of the differing RBE's (relative biological efficiency) of the various 
nuclear species. Finally, the theory does not yet include contributions from higher- 
order terms (above tertiary), which may be significant at greater absorber depths. 


CONCLUDING REMARKS 

Transport coefficients for heavy-ion transport in polyethylene and Du Pont 
Kapton have been calculated and utilized in an energy-dependent theory to predict 
depth-dose results for galactic, cosmic-ray, neon nuclei impinging upon a homogeneous 
shield composed of these organic composite materials. The importance of secondary 
and subsequent-generation ions produced by the nuclear fragmentation of the incident 
particles within the absorber was clearly demonstrated. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
December 12, 1983 
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SYMBOLS 


A j 

D ( x) 
E 

E x 

E 0 

j 

k 


ra jk (E) 


°m ( e) 

p m< e > 

Rj(E) 

rT 1 [R. (E) 1 
1 D 

Sj( E ) 

j (E) 


X 


Pi 

Oj (E) 

°ij (E) 

$ j (x,E) 


atomic mass of type-j ion, amu 
mass number, dimensionless 

energy absorbed per unit mass at x, MeV/g 

ion kinetic energy, MeV/amu 

residual energy, MeV/amu 

incident beam energy, MeV/amu 

type-j ion 

type-k ion 

multiplicity of type-j ions produced by collisions of type-k ions of 
energy E 

average extinction coefficient for type-M ions, cm ” 1 

nuclear attenuation factor for type-M ion at energy E, dimensionless 
continuous slowing-down range of type-j ion of energy E, cm 
inverse function of Rj(E) 

total stopping power or linear energy transfer (LET) due to interaction 
of type-j ion with orbital electrons of transport medium, MeV/cm 

stopping power or linear energy transfer (LET) per nucleon due to 

interaction of type-j ion with orbital electrons of transport medium, 
MeV/amu-cm 

o 

one-dimensional position vector, g/cm 

atomic number of type-j ion 

range scale parameter for type-j ion 

-3 

number density of ith constituent of absorber, cm 

macroscopic absorption cross section for type-j ion of energy E, cm 1 

microscopic absorption cross section for type-j particle of energy E 
colliding with type-i particle in absorber, cm Z 

differential flux of type-j ions at x with energy E, 

2 —1 
(cm -sec-MeV/amu) 



Subscripts : 


M type of ions in monoenergetic beam 

P projectile 

S Silberberg-Tsao formalism 

Superscripts : 

( 0 ) , ( 1 ) , ( 2 ) terms in series approximation to equation (3) 

Primes indicate a variable of summation or integration. 
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TABLE 1.- CONSTITUENT NUMBER DENSITIES FOR KAPTON 
AND POLYETHYLENE 


Constituent 

Number density, 
1 0 22 atoms/cm 3 , 

p i ' 
for - 

Kapton 
(1*42 g/cm 3 ) 

Polyethylene 
(0.91 g/cm 3 ) 

Hydrogen 

2.239 


7.826 

Carbon 

4.925 


3.913 

Nitrogen 

.448 



Oxygen 

1 .119 




TABLE 2.- MACROSCOPIC NUCLEAR-ABSORPTION CROSS SECTIONS IN KAPTON 


MeV /amu 


Macroscopic 
CTj (E) 

nuclear-absorption cross section, 
, 10” 4 /cm, in Kapton for - 


’h 

4 He 

7 Li 

9 Be 

12 c 

16 0 

20 Ne 

27 A1 

40 Ar 

56 Fe 

64 Cu 

25 

323 

526 

789 

862 

861 

1012 

1084 

1118 

1455 

1630 

1690 

50 

249 

430 

645 

713 

727 

858 

924 

971 

1260 

1430 

1490 

75 

215 

387 

582 

648 

667 

790 

853 

905 

1175 

1 340 

1400 

100 

191 

355 

534 

597 

621 

738 

799 

855 

1109 

1273 

1334 

200 

158 

308 

466 

526 

556 

663 

722 

782 

1016 

1176 

1237 

300 

149 

296 

447 

506 

538 

643 

701 

762 

990 

1150 

1210 

400 

152 

300 

452 

511 

543 

648 

707 

767 

997 

1156 

1217 

600 

175 

328 

492 

554 

583 

694 

754 

813 

1052 

1215 

1276 

1000 

189 

346 

517 

581 

609 

723 

784 

842 

1089 

1253 

1314 

2000 

191 

353 

525 

589 

618 

732 

794 

852 

1100 

1266 

1 330 

4000 

188 

350 

520 

584 

614 

728 

789 

848 

1095 

1260 

1 322 

10000 

185 

347 

515 

578 

609 

722 

783 

843 

1088 

1250 

1315 

























TABLE 3.- MACROSCOPIC NUCLEAR-ABSORPTION CROSS SECTIONS IN POLYETHYLENE 


[p = 0.91 g/cm 3 ] 


E, 

MeV/amu 

Macroscopic nuclear-absorption cross section, 
Cj(E), 10“ 4 /cm, in polyethylene for - 


'H 

4 He 

7 Li 

9 Be 

12 c 

16 0 

to 

o 

a 

(D 

27 A1 

O 

> 

n 

56 Fe 

64 Cu 

25 

256 

428 

717 

795 

786 

951 

1038 

1078 

1471 

1667 

1733 

50 

179 

342 

564 

635 

648 

787 

861 

919 

1242 

1433 

1500 

75 

145 

304 

499 

565 

585 

713 

783 

848 

1143 

1330 

1402 

100 

127 

275 

448 

511 

536 

655 

723 

792 

1067 

1250 

1323 

200 

105 

235 

376 

435 

468 

572 

636 

710 

956 

1 1 34 

1 206 

300 

100 

224 

357 

414 

448 

550 

613 

686 

922 

1100 

1172 

400 

104 

227 

361 

418 

454 

556 

619 

692 

930 

1110 

1180 

600 

125 

252 

400 

460 

496 

606 

671 

743 

994 

1177 

1248 

1000 

1 39 

268 

425 

488 

523 

638 

705 

775 

1036 

1 220 

1 292 

2000 

138 

274 

433 

497 

531 

646 

715 

786 

1050 

1236 

1 308 

4000 

134 j 

| 272 

427 

490 

527 

641 

708 

781 

1041 

1 228 

1 300 

10000 

132 

: 270 

421 

484 

522 

634 

_ 702 

774 

1032 

1220 

1290 


TABLE 4.- AVERAGE EXTINCTION COEFFICIENTS FOR IONS IN KAPTON 


E, 



Average 

10“ 

extinction coefficient, 
4 /cm, in Kapton for - 

V 



MeV/amu 

’h 

4 He 

7 Li 

9 Be 

12 c 

16 0 

20,_ 

Ne 

30 Si 

o 

> 

ri 

50y 

60 Ni 

25 

323 

526 

789 

862 

860 

1011 

1083 

1167 

1452 

1572 

1659 

50 

290 

483 

725 

795 

801 

944 

1013 

1102 

1371 

1490 

1579 

75 

259 

443 

665 

734 

746 

880 

947 

1040 

1 291 

1409 

1497 

100 

236 

414 

622 

689 

704 

833 

898 

993 

1232 

1347 

1435 

200 

189 

351 

528 

591 

615 

731 

793 

892 

1103 

1213 

1 300 

300 

171 

325 

491 

552 

580 

691 

750 

851 

1051 

1160 

! 1246 

400 

163 

315 

475 

535 

565 

673 

733 

833 

1029 

1 1 37 

1222 

600 

163 

314 

473 

534 

564 

672 

731 

832 

1026 

1134 

1220 

1000 

174 

327 

492 

553 

582 

693 

753 

854 

1052 

1 161 

1247 

2000 

184 

342 

510 

574 

602 

715 

776 

877 

1079 

1189 

| 1276 

4000 

187 

347 

517 

581 

610 

724 

785 

886 

1090 

1 200 

; 1287 

10000 

187 

348 

517 

581 

611 

724 

785 

887 

1091 

1200 

| 1287 





TABLE 5 .- AVERAGE EXTINCTION COEFFICIENTS FOR IONS IN POLYETHYLENE 


[p = 0.91 g/cm**] 





Average 

extinction 

coefficient. 

V 



E, 



10" /cm, in 

polyethylene for 




MeV/amu 

1 H 

4 He 

7 Li 

9 Be 

12 c 

16 0 

<D 

O 

CM 

30_. 

Si 

O 

> 

50 

V 

60 . 
Ni 

25 

256 

428 

717 

795 

786 

950 

1037 

1 1 35 

1468 

1603 

1700 

50 

221 

389 

648 

723 

725 

878 

959 

1064 

1372 

1506 

1605 

75 

189 

354 

586 

657 

667 

810 

886 

996 

1 279 

1410 

151 1 

100 

167 

328 

540 

609 

623 

759 

832 

945 

1209 

1338 

1439 

200 

128 

272 

442 

504 

530 

647 

715 

832 

1058 

1180 

1281 

300 

115 

250 

403 

463 

492 

602 

668 

786 

996 

1116 

1215 

400 

110 

240 

386 

445 

476 

583 

648 

765 

969 

1088 

1187 

600 

112 

240 

383 

442 

476 

582 

646 

764 

965 

1085 

1183 

1000 

123 

252 

400 

461 

495 

605 

671 

788 

994 

1115 

1 214 

2000 

133 

264 

419 

481 

516 

629 

696 

813 

1026 

1148 

1247 

4000 

135 

269 

425 

488 

523 

637 

705 

823 

1037 

1160 

1260 

10000 

133 

270 

424 

487 

524 

637 

704 

824 

i 1036 

1 

1160 I 

1260 


TABLE 6.- 20 Ne FRAGMENTATION PARAMETERS IN KAPTON AND 
TOTAL FRAGMENT CHARGE Z 


20 . 


10 


Ne fragmentation parameters in Kapton for 
MeV/amu, of - 




31 .6 


100 


316 


1000 


3160 


0 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 


1 .312 
3.1 39 
.971 
.036 
.043 
.003 
.005 
.020 
.193 
.302 
.008 
.043 

10.4 


0.968 

2.170 

.713 

.024 

.028 

.013 

.026 

.048 

.225 

.268 

.159 

.009 

10.2 


0.961 
1 .974 
.587 
.032 
.037 
.035 
.080 
.090 
.289 
.217 
.122 
.002 

10.2 


1 .287 
2.375 
.524 
.047 
.051 
.061 
.1 21 
.106 
.267 
.195 
.077 
.000 

10.2 


2.158 
3.177 
.651 
.083 
.080 
.074 
.1 1 1 
.083 
.183 
.158 
.060 
.000 

10.2 


3.420 
3.81 2 
.600 
.095 
.087 
.073 
.100 
.073 
.158 
.121 
.053 
.000 

10.0 




































20 . 


Ne fragmentation parameters in polyethylene for 
MeV/amu, of - 


E, 



i° 

31 .6 

100 

316 

1000 

3160 

0 

1 .313 

0.968 

0.962 

1 .288 

2.162 

3.427 

1 

3.141 

2.171 ; 

1 .975 

2.378 

3.182 

3.819 

2 

.971 

.713 

.588 

.525 

.653 

.601 

3 

.035 

.023 

.031 

.046 

.081 

.093 

4 

.042 

.026 

.036 

.050 

.078 

.086 

5 

.003 

.013 

.035 

.060 

.074 

.072 

6 

.005 

.026 

.080 

.1 21 

.111 

.100 

7 

.020 

.048 

.090 

.106 

.083 

.073 

8 

.193 

.225 

.289 

.267 

.183 

.159 


.303 

.268 

.217 

.195 

.158 

.121 


.008 

.159 

.122 

.077 

.060 

.053 

mm 

.043 

.009 

.002 

.000 

.000 

.000 

z 

F 

10.4 

10.2 

10.2 

i 

10.2 

10.2 

10.0 


















Particle energy, MeV/arcu 


Figure 1*- Calculated stopping powers in 
polyethylene for typical cosmic-ray 
ions as a function of kinetic energy. 



Figure 2.- Calculated stopping powers in 
Kapton for typical cosmic-ray ions as 
a function of kinetic energy. 



Relative dose 



Depth in absorber, cm 


Figure 3.- Relative dose as a function of 
depth in absorber for 350-MeV/amu 
^°Ne incident on Kapton. 



Figure 4.- Relative dose as a function of 
depth in absorber for 670-MeV/amu ^Ne 
incident on Kapton. 
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Depth in absorber, cm 


Figure 5.- Relative dose as a function of 
depth in absorber for 2000-MeV/amu z N 
incident on Kapton. 



o 5 10 15 20 

Depth in absorber, cm 


Figure 6*- Relative dose as a function of 
depth in absorber for 350-MeV/amu 
2( ^Ne incident on polyethylene. 



Relative dose ^ Relative dose 



igure 7.- Relative dose as a function of 
depth in absorber for 670-MeV/amu 
^Ne incident on polyethylene. 



Figure 8.- Relative dose as a function of 
depth in absorber for 2000-MeV/amu 
20 Ne incident on polyethylene. 
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